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ABSTRACT 


Current  active  imaging  algorithms  for  moving  targets  suffer  from  issues  of 
incorrect  positions  (spatial)  and  streaking  artifacts  (temporal).  Using  the  Cheney/Borden 
procedure,  we  investigated  combining  the  spatial,  temporal,  and  spectral  aspects  of  real 
and  synthetic  aperture  radar  images.  We  code  the  Cheney/Borden  algorithm  to  include 
the  target  velocity,  include  an  appropriate  threshold,  and  illustrate  how  multistatic  radar 
can  determine  a  target’s  location  in  phase  space.  By  running  simulations  on  single  and 
multiple  moving  targets,  we  showed  that  an  iteration  of  velocity  and  position  choices  for 
targets  enhanced  the  correlation  map  for  multistatic  radar  systems. 


V 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


VI 


TABLE  OF  CONTENTS 


L  INTRODUCTION . 1 

A.  INTRODUCTION . 1 

B.  MAXWELL’S  EQUATIONS . 1 

C.  ELECTROMAGNETIC/RADAR  THEORY . 2 

1.  Electromagnetic  Theory . 2 

2,  Radar  Theory . 3 

D.  RADAR  RANGE-PROFILES . 3 

E.  MULTIPLE-LOOKS . 4 

F.  FILTERED  BACKPROJECTION . 5 

G.  SYNTHETIC  APERTURE  IMAGING . 6 

H.  MOVING  TARGETS . 7 

II.  IMAGE  ALGORITHM  PROGRESSION . 9 

A,  THE  CHENEY/BORDEN  ALGORITHM . 9 

1,  Scattered  Field  for  Moving  Targets . 9 

2,  Imaging  via  a  Filtered  Adjoint . 10 

3,  Image  Analysis . 11 

B,  TAN  LU  PIN’S  WORK . 11 

C,  TEOBENGKOON  WILLIAM’S  WORK . 12 

III.  DATA  ANALYSIS . 13 

A.  STRUCTURE . 13 

B,  CODING . 13 

1,  Cheney/Borden  Code . 14 

2,  Thresholding . 14 

3,  Normalization . 15 

4,  Display . 15 

IV.  FINDINGS  AND  FUTURE  WORK . 17 

A,  FINDINGS . 18 

B,  FUTURE  WORK . 22 

V.  SUMMARY  AND  CONCLUSION . 25 

A,  SUMMARY . 25 

B.  CONCLUSION . 25 

APPENDIX . 27 

LIST  OF  REFERENCES . 31 

INITIAL  DISTRIBUTION  LIST . 33 


vii 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


LIST  OF  FIGURES 


Figure  1.  Electromagnetic  spectrum  (From  [1]) . 1 

Figure  2.  Basic  radar  system  (From  [2]) . 3 

Figure  3.  Single  pulse  ambiguity  for  targets  of  the  same  range  (From  [2]) . 4 

Figure  4.  Multiple  Pulse  Radar  (From  [2]) . 4 

Figure  5.  Filtered  backprojection  shown  using  single  to  multiple  view  angles 

(From  [2]) . 5 

Figure  6.  SAR  and  ISAR  schemes  for  imaging  targets  (From  [2]) . 6 

Figure  7.  In  spotlight  SAR  only  one  area  is  imaged  (From  [2]) . 6 

Figure  8.  Stripmap  SAR  acquires  a  sequence  of  radar  returns  along  a  path  (From 

[2]) . 7 

Figure  9.  Display  progressions  of  target  located  (-50,0)m  with  one  transmitter 
located  (-500,-500)m  and  four  receiver  located  at  (-500,-500)m, 
(-500,500)m,  (500,500)m,  (500,-500)m.  (a)  raw  data  (b)  thresholded  data 

(c)  normalized  data  (d)stem/quiver  plot  of  data,  (x-y  axis  are  in  meters) . 13 

Figure  10.  Max  amplitude  of  velocity  profiles  at  15  degree  cross-ranges . 14 

Figure  11.  (a)One  transmitter  (O)  and  1  receiver  (X),  (b)  1  transmitter  and  4 
receivers,  (c)  1  transmitter  and  11  receivers,  (d)  1  transmitter  and  21 

receivers . 17 

Figure  12.  (a)  imagesc  and  (b)  surf  display  of  single  target,  transmitter  and  receiver 

(x-y  axis  are  in  meters) . 18 

Figure  13.  (a)  imagesc  and  (b)  surf  display  of  1  transmiter  and  4  receivers  arranged  as 

in  Figure  11b  (x-y  axis  are  in  meters) . 18 

Figure  14.  (a)  imagesc  and  (b)  surf  display  of  Itransmiter  and  4  receivers  arranged  as 

in  Figure  11b  (x-y  axis  are  in  meters) . 19 

Figure  15.  (a)  imagesc  and  (b)  surf  display  Itransmiter  and  11  receivers  arranged  as 

in  Figure  lie  (x-y  axis  are  in  meters) . 19 

Figure  16.  (a)  imagesc  and  (b)  surf  normalized  display  Itransmiter  and  21  receivers 

arranged  as  in  Figure  lid  (x-y  axis  are  in  meters) . 20 

Figure  17.  (a)  &  (b)  surf  plots  of  the  same  targets  where  (a)  has  targets  of  equal 

velocity,  (b)  has  one  target  velocity  significantly  greater  than  the  other  (x- 

y  axis  are  in  meters) . 20 

Figure  18.  Multiple  targets  with  a  threshold  applied,  at  significantly  different 

velocities  (x-y  axis  are  in  meters) . 21 

Figure  19.  Plot  of  five  targets  using  Matlab  stem  and  quiver  located  (0,0)m, 

(50,50)m,  (50,-50)m,  (-50,50)m,  (-50,-50)m  with  different  velocities  (x-y 
axis  are  in  meters) . 22 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


X 


I.  INTRODUCTION 


A,  INTRODUCTION 

This  first  chapter  is  an  overview  of  the  physicals  foundations  used  in  the 
simulations  leading  to  synthetic  aperture  radar.  This  overview  is  presented  as  a 
foundation  for  Chapter  II,  which  reviews  previous  work. 


B.  MAXWELL’S  EQUATIONS 
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Figure  I.  Electromagnetic  spectrum  (From  [1]) 


Observations  of  the  electromagnetic  spectrum  (Figure  1)  and  the  associated 
Forentz  force  law  provided  physicists  with  the  building  blocks  required  to  explain  the 
electromagnetic  spectrum  in  terms  of  two  fields.  The  fields  are  the  Electric  (E)  and 
Magnetic  (B)  fields  and  they  determine  the  forces  felt  by  any  charged  object.  In  free 
space,  monochromatic  versions  of  each  field  travel  in  a  sinusoidal  pattern;  they  are 
perpendicular  to  each  other  and  are  in  phase  with  each  other.  These  observations  led 
physicists  to  develop  the  building  blocks  for  the  Maxwell  equations. 

James  Clerk  Maxwell  was  known  for  many  things,  but  his  greatest 
accomplishment  was  completing  and  consolidating  electro-magnetic  theory.  Maxwell 
realized  that  Gauss’  law  for  Electricity, 

VD  =  p  (1.1) 

where  D  =  sE  ,  and  p  is  the  charge  density,  the  equivalent  law  for  magnetism 

V-5  =  0  (1.2) 


1 


Faraday’s  law, 


VxE  = - 

dt 


and  Ampere’s  law  with  Maxwell’s  eorreetion 


VxH  =  J  +  — 
dt 


(1.3) 


(1.4) 


where  H  =  Bj /u,  and  by  the  inelusion  of  the  displacement  current  J  could  be  used  to 
define  the  four  fundamental  laws  of  Electromagnetic  theory,  and  the  theory  of  light. 


C.  ELECTROMAGNETIC/RADAR  THEORY 

The  principles  of  electromagnetic  and  radar  theory  rely  on  the  Maxwell's 
equations.  Radar  theory  is  a  practical  expansion  of  the  fundamental  theory  of 
electrodynamics  and  is  implemented  in  many  radar  systems  in  use  today. 


1.  Electromagnetic  Theory 


A  key  building  block  to  electromagnetic  theory  is  the  electromagnetic  wave 
equation.  By  simple  mathematical  manipulation  Maxwell’s  equations  become 
inhomogeneous  wave  equations. 


'  ^^2  A*  '  "  J  • 

Ot 

In  free  space  they  become  the  homogeneous  wave  equation, 

d^B  „ 


The  same 


dt 

manipulations  applying  to  the  electric  field  give: 

d^E  dJ  Vp 
V  E-pS - 2r  =  p - + - 

^  dt  S 


These  wave 


V  r.  —  Ub - ^  —  u - 1 - 

de  dt  a 

equations  provide  the  basic  foundation  for  numerous 


(1.5) 

(1.6) 

(1.7) 
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2. 


Radar  Theory 


target 


Figure  2.  Basic  radar  system  (From  [2]) 


Figure  2  illustrates  electric  fields  and  currents  induced  on  an  antenna  by  a 
transmitter.  Those  fields  produce  an  electromagnetic  pulse  that  propagates  through  the 
medium  and  eventually  interacts  with  a  target.  The  electromagnetic  pulse  induces  a 
current  in  the  target,  which  in  turn  creates  an  electromagnetic  pulse  in  response.  This 
“scattered”  electromagnetic  pulse  propagates  back  through  the  medium  to  the  antenna 
where  it  induces  a  field.  Using  theories  developed  from  the  wave  equations,  the  range 
can  be  determined  from  the  signal  resulting  from  the  return  field  induced  on  the  antenna. 

D.  RADAR  RANGE-PROFILES 

Current  radars  assume  targets  are  effectively  stationary.  This  is  a  good 
approximation  since  the  duration  of  the  pulses  are  relatively  small  and  the  speed  of  the 
targets  are  a  small  fraction  of  the  speed  of  light.  Range-profiles  are  one-dimensional  high 
range  resolution  images  of  the  target  using  the  assumption  targets  are  stationary.  It  is 
difficult  to  use  range-profiles  for  target  identification  since  all  the  scatterers  located  at  the 
same  distance  from  the  receiver  have  the  same  time  delay.  A  single  range-profile  will  not 
be  able  to  distinguish  a  cross-range  structure. 
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E.  MULTIPLE-LOOKS 
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Figure  3.  Single  pulse  ambiguity  for  targets  of  the  same  range  (From  [2]) 

Figure  3  illustrates  how  a  single-range  profile  is  not  able  to  distinguish  a  eross- 
range  structure  of  multiple  targets  at  equal  ranges. 


Multiple  pulse  “imaging 


o 


o 


o 


Figure  4.  Multiple  Pulse  Radar  (From  [2]) 


A  multiple  pulse  radar  exploits  multiple  pulse  data  collected  from  multiple  aspect 
angles  to  distinguish  multiple  targets  that  are  at  the  same  range.  By  applying  the  simple 
concept  of  triangulation,  a  multiple  aspect  radar  is  able  to  determine  the  range  and  cross¬ 
range  of  multiple  targets.  The  multiple-look  concept  is  a  crude  illustration  of  “filtered 
backprojection.” 
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F. 


FILTERED  BACKPROJECTION 


8  views  60  views 

Figure  5.  Filtered  baekprojeetion  shown  using  single  to  multiple  view  angles 

(From  [2]) 


Filtered  baekprojeetion  takes  multiple  aspect  radar  range  profiles,  aligns  them  to  a 
common  origin,  and  sums  the  aligned  data  to  produce  an  image  of  the  scanned  area.  As 
the  number  of  looks  at  multiple  angles,  in  Figure  5,  increases  from  1  to  8,  a  discemable 
image  starts  to  appear  from  the  alignment  and  summation  of  the  data.  With  60  views  at 
different  angles,  the  filtered  baekprojeetion  of  the  data  displays  a  clear  image.  This 
provides  a  method  to  display  multiple  looks  of  data  using  the  fdtered  baekprojeetion 
method. 
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G.  SYNTHETIC  APERTURE  IMAGING 


SAR 


Figure  6.  SAR  and  ISAR  schemes  for  imaging  targets  (From  [2]) 


Synthetic  aperture  imaging  can  be  accomplished  using  a  stationary  antenna  and 
rotating  target  or  a  stationary  target  with  a  moving  antenna.  Synthetic  aperture  radar 
(SAR)  uses  the  concept  of  a  moving  antenna  that  illuminates  a  stationary  target  with  a 
series  of  pulses  from  the  moving  antenna  (Figure  6).  Inverse  synthetic  aperture  radar 
(ISAR)  assumes  a  stationary  radar  radiating  a  moving  target  (Figure  6). 

SAR  imaging  can  be  done  either  by  tracking  a  single  target  of  interest,  known  as 
spotlighting  or  by  stripmapping. 


Figure  7.  In  spotlight  SAR  only  one  area  is  imaged  (From  [2]) 
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In  spotlight  SAR  a  platform  illuminates  a  specific  area  to  get  multiple  images  of 
the  area  from  multiple  aspect  angles  on  the  area.  Spotlight  SAR  imaging  is  illustrated  in 
Figure  7. 


Figure  8.  Stripmap  SAR  acquires  a  sequence  of  radar  returns  along  a  path 

(From  [2]) 

Stripmap  SAR  requires  a  platform  to  sequentially  illuminate  a  moving  patch  and 
then  combines  the  measured  returns  to  synthesize  a  larger  image  than  typically  available 
from  spotlight  SAR. 

H,  MOVING  TARGETS 

When  a  target  moves,  it  rarely  travels  in  a  predictable  manner.  Unpredictable 
targets  can  cause  problems  with  synthetic  aperture  radars  because  current  SAR 
algorithms  assume  a  stationary  target.  Moving  targets  add  doppler  frequency  and  phase 
shift  artifacts.  The  most  notable  of  problems  is  image  artifacts.  Image  artifacts  come 
from  many  sources:  tank  cannon  barrels,  jet  engines,  or  (the  most  notable)  the  “train  off 
the  track”  effect.  The  Cheney/Borden  algorithm  [3]  does  not  assume  a  stationary  target 
and,  hopefully,  images  a  moving  train  that  is  on  the  tracks. 
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II.  IMAGE  ALGORITHM  PROGRESSION 


Chapter  II  will  survey  work  previously  eompleted,  which  will  explain  the 
chronology  of  developments  of  the  Cheney/Borden  algorithm. 

A,  THE  CHENEY/BORDEN  ALGORITHM 

The  Cheney/Borden  Algorithm  consists  of  a  scattering  model,  imaging  via  a 
filtered  adjoint,  and  image  analysis. 


1.  Scattered  Field  for  Moving  Targets 

Using  the  Born  Approximation  [4]  for  non-moving  targets  gives  the  scattered 

field  as 

'^scan  {x,t)  =  jjj  V ,  (2.1) 

D 


where  is  the  time-domain  Green  function  and  incident 

field.  It  can  be  shown[4]  that  when  the  incident  field  is  formed  using  a  time  domain 
signal  transmitted  from  position  y,  that  the  scattering  model  becomes: 


- 

■  Z  -X 

'/c)  Sy(t'  +  Ty 

x  - 

-y 

A) 

An; : 

-  — ^ 

^  -X 

An; 

X 

-y 

pi^x')d^x'dt' 


(2.2) 


where  >o(x')  is  the  spatial  scatterer  density. 


Tan  Lu  Pin  and  Teo  Beng  Koon  William  [4], [5], [6]  have  shown  that  a  slowly 
moving  target  illuminated  by  a  narrow-band  waveform  gives  an  approximate  form  for  the 
scattered  field  of  a  moving  target  of  the  form: 

ff  f  r-  -  -  it 

¥scaniy^^d')  ^  x2 1_|  I  _|  JJ  {y-z)-[x-v{t'  +  z-  x/c)J} 

(4;r)  p||t|  (2.4) 

+  [^y-zyxlc^p^[^x)d^vd^x  . 

In  this  expression,  we  can  ignore  the  constants  in  front  of  the  integrand  since  they 
are  not  needed  for  imaging.  Now,  assuming  the  field  is  from  an  unknown  point  scatterer 
located 
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at  position  p  and  velocity  u  ,  then  Equation  2.3  reduces  to  [4] 

xl^,[t'  +  {y  +  z)-plc)  . 

The  next  step  is  to  use  the  scattered  and  incident  fields  to  create  an  imaging 
function. 

2.  Imaging  via  a  Filtered  Adjoint 

Using  the  scattered  field  (Equation  2.2),  a  cross  correlation  gives  an  image  of  the 
form  [4] 

I{p,u)  =  \\\ [y, z, [y, z, t') dt'd’^yd’'z ,  (2.5) 

where  m,  n  =  1, 2, 3  depending  on  the  placement  of  the  transmitters  and  receivers. 

Because  there  is  a  chance  of  image  artifacts  for  moving  targets,  we  need  to  add  a 
filtering  function  Q[co,t',p,u,y,z^  so  that  [4] 

/ (^, M )  =  JJJ [y,  z, t')//*  [y, z,  t') Q[o), t' , p, u, y,z)dt'd”'yd’'z  .  (2.6) 

Expanding  Equation  2.6,  the  imaging  equation  becomes  [4] 

I(P,u)  =  -fjf 

jjj  (2.7) 

X  +  {y  +  z)-  pIc)ii/^^^„  [y, z, t') dt'd"yd"z 

Equation  2.7  forms  the  basis  of  an  imaging  algorithm. 
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3,  Image  Analysis 

Choosing  values  of  Q[co,t',p,u,y,z^  that  remove  the  |z|  and  \y\  coordinates 
from  the  amplitude  factor  in  Equation  2.3; 


- ^  (4;z-)|z||y|  _ 

p,u,y,z)  =  ^ — - '-:i[p,u,y,z), 

(o.. 


Q{co,t',p 

where  J(j?,M,j),z)  is  a  Jacobian  introduced  during  the  analysis.  Equation  2.7  becomes 
[4] 

/( =  J|  A(  j),z) p{x^d^vd^x ,  (2.8) 

where  p{x)  represents  a  perfect/ideal  image  function,  and  [4] 


K(^p,u,y,z)  =  ||exp|zA:^,  (j)  +  z)-[/>-x-M(z-/>)/c  +  v(z-x)/c]| 

X  j 5*  (t'  +  (y  +  z)  •  p/c)sy  (t'  +  ( j)  +  z)  •  x/c)  ^2.9) 

x](^p,u,y,z^d'"yd"z 

is  the  point  spread  function  describing  the  image. 

Making  a  change  of  variables,  Equation  2.9  becomes  [4], [5], [6] 

A(^,M,y,z)  =  JJexp{zl^,(j)  +  z)-[M(z-^)/c-v(z-x)/c]} 

xexp{zl^(y  +  z)-[(^-x)  +  i(M-v)(y  +  z)-(^  +  x)]/c}  (2.10) 

x2r(^>.(E  +  ^)-(w-v),(j)  +  z)-(^-x)/c)j(^,M,y,z)J'”yJ"z  , 
where  x  is  the  radar  ambiguity  function 


x(v,t)=  j  s*(t  +  jT)s^,(t-jr)e-‘''‘dt. 


B,  TAN  LU  PIN’S  WORK 

Tan  Eu  Pin  (Tan)  conducted  simulations  on  the  Cheney/Borden  algorithm 
(Equation  2.8).  Tan’s  work  showed  that  the  imaging  point  spread  function  is  well 
behaved,  localizes  the  target  in  phase  space  and  is  translation  invariant  [6].  Phase  space, 
in  this  context,  is  a  space  where  position  and  velocity  of  an  image  point  are  represented, 
with  correspond  to  one  unique  point  in  velocity  and  position  space. 
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Localizing  a  target  in  phase  space  demonstrated  that  a  moving  target  eould  be 
imaged  in  a  multistatie  environment  in  whieh  target  motion  is  not  ignored.  Tan  examined 
a  single  point  seatterer  in  a  2-D  plane. 

C.  TEOBENGKOON  WILLIAM’S  WORK 

Teo  Beng  Koon  William  (William)  used  the  Cheney/Borden  algorithm  (Equation 
2.8)  to  show  that  the  seattering  model,  radar  model,  and  imaging  model  eould  be 
eombined.  William  ran  simulations  on  each  part  and  proved  that  for  a  single  point 
seatterer  a  stationary  target’s  loeation  eould  be  determined  using  a  multistatie  array  [6]. 
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III.  DATA  ANALYSIS 


A.  STRUCTURE 

The  current  investigation  was  performed  by  Chee  Young  Ng,  Christopher  Carroll, 
and  the  author  and  started  where  Teo  Beng  Koon  William  concluded.  Chapter  III 
presents  the  chronology  of  work  conducted  utilizing  the  Cheney/Borden  algorithm 
showing  that  for  a  single  point  scatterer  a  multistatic  array  can  image  a  moving  target. 

B,  CODING 
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Figure  9.  Display  progressions  of  target  located  (-50,0)m  with  one  transmitter  located  (- 
500,-500)m  and  four  receiver  located  at  (-500,-500)m,  (-500,500)m,  (500,500)m, 
(500,-500)m.  (a)  raw  data  (b)  thresholded  data  (c)  normalized  data  (d)stem/quiver  plot  of 

data,  (x-y  axis  are  in  meters) 


A  Matlab  code  utilizing  the  Cheney/Borden  algorithm  appears  in  the  Appendix. 
This  code  provided  the  basis  for  a  series  of  simulation  runs.  The  code  has  four  parts. 
Including  velocity  and  verifying  imaging  of  a  moving  target  was  the  driving  force  behind 
coding.  Thresholding  was  needed  to  eliminate  the  image  artifacts  and  thereby  prevent 
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false  targets.  Normalizing  was  used  to  display  all  velocities  with  the  same  peak  value. 
Matlab  was  used  to  run  simulations  and  display  images  created  from  the  Cheney/Borden 
Algorithm. 

1,  Cheney/Borden  Code 

Up  until  now,  the  algorithm  was  run  without  considering  target  velocity  and 
proved  only  that  a  point  scatter  could  be  located  in  space.  Our  first  goal  included 
velocity  into  a  2D  code. 

The  Cheney/Borden  algorithm  produces  an  image  by  comparing  the  scattered 
field  of  a  known  target  to  the  measured  scattered  field.  The  trial  target’s  position  and 
velocity  must  be  specified  with  respect  to  each  transmitter  and  receiver  location.  Cross¬ 
correlation  forms  the  image.  In  our  testing  the  Matlab  code  imagesc  displays  the 
correlation  map  (which  constitutes  the  “image”)  (Figure  9a). 
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Figure  10.  Max  amplitude  of  velocity  profiles  at  15  degree  cross-ranges 


The  peaks  of  the  correlation  map  for  stationary  to  slow  moving  targets  were 
stronger  than  those  of  faster  moving  targets  (Figure  10),  causing  the  fast  movers  to  be  lost 
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in  the  artifacts.  The  imagesc  function  in  Matlab  showed  that  the  data  needed  to  be 
thresholded  to  decrease  the  number  of  false  returns  (Figure  9a).  The  space,  or  lack  there¬ 
of,  around  the  target  created  the  false  artifacts. 

Computing  the  minimum  amplitude  of  velocity  for  a  known  target  provided  the 
threshold  level.  A  conditional  in  the  code  set  any  computed  point  that  fell  below 
minimum  amplitude  for  the  given  velocity  to  zero.  By  doing  this,  simulations  with  no 
targets  moving  at  the  checked  velocities  were  set  to  zero.  Figure  9b  shows  the  data  with 
a  threshold  partition.  It  is  clear  that  the  artifacts  have  been  eliminated. 

3.  Normalization 

After  thresholding  the  correlation  map  image,  amplitude  normalization  allows  the 
imaging  of  multiple  targets  to  be  shown  with  the  same  peak  level  (Figure  9c).  A  partition 
in  code  allowed  the  data  to  be  normalized  for  each  velocity.  Matlab  code  imagesc  and 
surf  displayed  the  data  (Figure  9c). 

4,  Display 

Finally,  the  imaged  data  needed  to  display  the  position,  velocity  and  direction  of 
the  target.  The  ease  of  programming  and  graphing  in  Matlab  made  image  display 
straightforward.  The  Matlab  subroutine  imagesc  displays  the  data  in  a  2-D  format,  which 
provides  the  targets  position.  Matlab  subroutine  surf  shows  the  intensity  of  the  return 
from  the  imaged  target  in  the  form  of  a  peak.  The  Matlab  plot  functions  stem  with  quiver 
were  the  most  intuitive  with  respect  to  a  user-friendly  display  (Figure  9d). 
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IV.  FINDINGS  AND  FUTURE  WORK 


X  1o"  (a)  X  1o"  (b) 


x(m)  x(m) 

Figure  11.  (a)One  transmitter  (O)  and  1  reeeiver  (X),  (b)  1  transmitter  and  4  reeeivers,  (e)  1 
transmitter  and  11  reeeivers,  (d)  1  transmitter  and  21  reeeivers. 

Simulations  were  run  using  both  monostatie  and  multistatie  configurations  (see 
Figure  11).  The  transmitter  and  receivers  were  stationary  for  the  simulations.  The  units 

in  the  x-y  plane  are  in  meters  and  the  velocities  in  component  form  ,  are  in  meters 

per  second.  The  imaged  area  on  the  x-y  plane  is  a  square  of  dimension  200m  x  200m 
centered  at  the  origin  (chosen  to  reduce  Matlab  computation  time).  A  single  pulse  is  of 
Q.2/US  duration  with  a  100/^5  period.  A  pulse  train  is  three  Q-l/us  pulses  with  a  100/^5 
period. 
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Figure  12.  (a)  imagesc  and  (b)  surf  display  of  single  target,  transmitter  and  reeeiver  (x-y  axis 

are  in  meters) 


Figure  12a,  illustrates  an  image  formed  from  the  sensor  eonfiguration  of  Figure 
10a.  The  target  was  plaeed  at  (-50,0)m  moving  along  the  x-axis  with  veloeity  (10,0)m/s. 
A  single  pulse  produees  the  eorrelation  map  of  Figure  12  for  a  single  transmitter  and 
reeeiver  this  gives  the  range  but  not  the  eross-range  loeation  of  the  target. 


The  next  imaging  step  added  more  reeeivers  for  a  multistatie  look  at  the  area.  A 
single  pulse  illuminated  a  target  loeated  at  (-50,0)m  with  veloeity  (10,0)m/s.  The 
transmitter  and  reeeiver  setup  was  as  in  Figure  11b.  Figure  13a  shows  an  image  of  the 
target  from  the  data,  but  there  are  image  artifaets  that  eonfound  image  interpretation. 
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Figure  14.  (a)  imagesc  and  (b)  surf  display  of  Itransmiter  and  4  reeeivers  arranged  as  in 

Figure  11b  (x-y  axis  are  in  meters) 


The  image  in  Figure  14  uses  the  same  set-up  for  the  target  and  reeeivers  as  in 
Figure  13.  Figure  14  uses  a  pulse  train  for  the  image  redueing  the  number  of  image 
artifaets. 
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Figure  15.  (a)  imagesc  and  (b)  surf  display  Itransmiter  and  1 1  reeeivers  arranged  as  in 

Figure  lie  (x-y  axis  are  in  meters) 
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Figure  15  shows  the  target  at  (-50,0)m  with  a  veloeity  of  (10,0)m/s  and  imaged 
using  a  pulse  train.  Figure  15  shows  that  by  inereasing  the  number  of  reeeivers  the  image 
quality  inereases.  Slower  target  speeds  produee  greater  correlation  peaks  (Figure  15b). 
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Figure  16.  (a)  imagesc  and  (b)  surf  normalized  display  Itransmiter  and  21  receivers 

arranged  as  in  Figure  lid  (x-y  axis  are  in  meters) 


Figure  16  reveals  that  increasing  the  number  of  receivers  looking  at  an  area 
produces  a  clearer  image  and  reduces  the  image  artifacts. 


Figure  17.  (a)  &  (b)  surf  plots  of  the  same  targets  where  (a)  has  targets  of  equal  velocity,  (b) 

has  one  target  velocity  significantly  greater  than  the  other 
(x-y  axis  are  in  meters) 
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Figure  17  incorporates  targets  with  different  velocities.  Figure  17a  has  targets 
located  at  (-50,0)m,  (50,0)m  with  a  velocities  of  (10,0)m/s,  (-10,0)m/s  respectfully,  and 
Figure  17b  has  targets  at  (-50,0)m,(50,0)m  with  velocities  (-10,0)m/s,(-100,0)m/s 
respectfully.  When  the  targets  move  at  different  velocities  but  near  the  same  speed 
(Figure  17a),  the  peaks  are  well  above  the  artifacts.  Figure  17b  shows  that  for 
significantly  different  speeds  greater  than  a  factor  of  10,  the  faster  target  is  lost  in  the 
artifacts.  It  was  observed  that  thresholding  would  be  required  to  see  multiple  targets 
moving  at  significantly  different  speeds.  Thresholding  would  also  be  needed  to  remove 
the  image  artifacts  from  a  correlation  map  with  no  target  moving  at  the  velocity  used  to 
compute  the  algorithm  (from  Chapter  111). 


Thresholding 


Figure  18.  Multiple  targets  with  a  threshold  applied,  at  significantly  different  velocities  (x-y 

axis  are  in  meters) 

Figure  18  displays  the  thresholded  version  of  the  data  used  in  Figure  17.  The 
thresholding  was  done  by  inserting  a  loop  into  Matlab  and  a  peak  conditional  that  found 
the  peak  of  the  signal  and  removed  the  artifacts  by  setting  a  limit  that  an  expected  peak 
must  exceed.  When  the  data  was  summed  and  normalized,  both  peaks  remained  resulting 
in  the  image  of  the  positions  for  each  target. 
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Figure  19.  Plot  of  five  targets  using  Matlab  stem  and  quiver  located  (0,0)m,  (50,50)m,  (50,- 
50)m,  (-50,50)m,  (-50,-50)m  with  different  velocities  (x-y  axis  are  in  meters) 

Figure  19  illustrates  that  the  algorithm  can  distinguish  multiple  targets  moving  at 
multiple  velocities  in  2-D  space.  There  were  five  targets  moving  with  different  velocities 
with  the  transmitter  and  receivers  setup  as  in  Figure  lOd  for  this  simulation.  The 
Cheney/Borden  algorithm  was  run  with  the  image  being  produced  by  the  Matlab  codes 
stem  and  quiver.  Stem  provided  the  location  of  the  target  while  quiver  produced  the 
velocity  magnitude  arrows.  Using  this  setup  a  user  would  be  able  to  see  where  the  target 
is  located  and  what  velocity  it  has. 

The  targets  for  Figure  19  are  point  targets  in  a  2-D  world.  For  actual  data  targets, 
the  data  would  be  a  grouping  of  points  with  arrows  from  each  point  pointing  in  the  same 
general  direction. 

B,  FUTURE  WORK 

Future  work  on  the  project  will  be  in  five  areas.  The  first  area  is  in  the  speed  at 
which  the  program  runs.  Matlab  is  a  user-friendly  interface  that  is  easy  to  use  but  is  it 
slow.  Each  simulation  run  took  on  average  40+  minuets  to  run.  Taking  the  code  and 
programming  it  into  C  as  a  .mex  file  to  be  called  by  Matlab  would  greatly  increase  the 


stem  &  quiver 
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speed  of  the  program.  A  faster  program  would  greatly  increase  the  number  of 
simulations  that  could  be  run  in  a  given  day. 

Being  able  to  run  more  simulations  would  allow  more  display  options  to  be 
considered.  Displaying  the  data  in  a  user-friendly  manner  is  essential  if  anyone  other 
than  someone  familiar  with  the  code  is  to  interpret  it. 

In  principle,  the  fast-slow  velocity  correlations  should  be  of  the  same  magnitude, 
independent  of  velocities.  We  suspect  that  poor  resolution  produces  poor  correlation  for 
larger  velocities.  Solving  the  correlation  issues  for  large  velocities  should  provide  a 
better  threshold  value. 

A  big  step  in  the  progress  of  the  project  is  to  apply  the  algorithm  to  real  data. 
Procuring  unclassified  data  with  which  to  run  simulations  is  vital.  Obtaining  such  data  is 
one  step  of  this  project,  while  another  step  would  be  in  assimilating  the  data  into 
something  that  could  be  run  with  the  algorithm. 

While  running  the  algorithm  an  appropriate  form  for  the  filtering  function 
Q[co,t',p,u,y,z^  is  needed.  For  the  current  project,  the  filtering  function  used  1  in 

Equation  2.7  for  expedience.  Depending  on  the  medium  and  time  of  travel,  the 
algorithm  will  need  to  be  run  multiple  times  to  determine  a  best  fit  filtering  function, 
bringing  up  the  subject  again  of  needing  a  faster  way  to  run  the  algorithm. 
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V.  SUMMARY  AND  CONCLUSION 


A.  SUMMARY 

In  this  thesis,  we  found  that  the  Cheney/Borden  algorithm  eould  image  a  moving 
target  and  reduee  the  image  artifaets  assoeiated  with  moving  target  images.  It  was  found 
that  eaeh  time  the  algorithm  was  run  for  a  speeified  veloeity,  the  targets  moving  at  that 
veloeity,  eould  stand  out  in  the  displayed  image.  For  a  speeified  range  of  veloeities,  the 
program  ean  be  restrieted  to  eonsider  only  targets  falling  within  that  veloeity  range  for 
imaging  purposes.  For  example,  for  images  on  the  ground  a  vehiele,  boat,  or  other  target 
on  the  ground  should  not  be  moving  faster  than  70  miles  per  hour.  The  program  ean 
restriet  target  images  to  objeets  between  0  to  70  miles  per  hour.  Constraining  the  veloeity 
spaee  will  allow  the  images  to  be  ereated  faster  allowing  for  faster  review  of  the  images. 

B,  CONCLUSION 

The  Cheney/Borden  algorithm  is  a  tool  that  ineorporates  the  spatial,  temporal  and 
speetral  aspeets  of  radar  returns  into  one  equation.  By  eombining  these  three  aspeets, 
multiple  radar  systems  may  no  longer  be  required  to  image  an  area.  Being  able  to  image 
an  objeet  while  simultaneously  determining  its  position,  heading,  and  veloeity  will 
revolutionize  the  radar  eommunity.  The  projeet  is  one  step  eloser  by  eonfirming  that  a 
multistatie  array  ean  image  an  area  and  determine  the  position  and  veloeity  of  multiple 
moving  targets. 
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APPENDIX 


The  Cheney/Borden  algorithm  was  coded  in  Matlab,  and  shows  how  normalizing, 
thresholding,  and  displaying  of  the  data  is  preformed  to  provide  an  image  of  the  type  in 
Chapters  IV  and  V. 


%Origin  of  the  image  scene  is  set  at  (0,0) 

%Target  Information 

N_tt=l;  %No  of  Targets 

tt=[-50  0;50  0];  %Target  X  -  Y  postion 

tt  vel=[l  -25; -75  -10];  %Target  velocity  in  X  -  Y  direction 

%Transmitter  Information 

N_Tx=l;  %No  of  Transmitter 

Tx  =[0  -10e3];  %Transmitter  X  position,  Y  position 
T  tx=0;  %Start  time  of  transmitted  pulse 

Tx  mag=  sqrt (Tx ( 1 , 1 ) ^2+Tx ( 1 , 2 ) ^2 ) ; 

Tx_hat  =  Tx/sqrt (Tx (1, 1) ^2+Tx (1, 2) ^2) ; 

%Receiver  Information 
N  Rx=21;  %number  of  receivers 

Rx=[-10e3  -10e3;-9e3  -10e3;-8e3  -10e3;-7e3  -10e3;-6e3  -10e3;-5e3  - 
10e3;-4e3  -10e3;-3e3  -10e3;-2e3  -10e3;-le3  -10e3;0  -10e3;le3  -10e3;2e3 
-10e3;3e3  -10e3;4e3  -10e3;5e3  -10e3;6e3  -10e3;7e3  -10e3;8e3  -10e3;9e3  - 
10e3;10e3  -10e3] ;  %Receiver  X  -  Y  position 
Rx  mag=  sqrt (Rx ( : , 1 ) . *Rx ( : , 1 ) +Rx ( : , 2 ) . *Rx ( : , 2 ) ) ; 

Rx_hat ( : , 1 ) =Rx ( : , 1 ) . /sqrt (Rx ( : , 1 ) . *Rx ( : , 1 ) +Rx ( : , 2 ) . *Rx ( : , 2 ) ) ; 

Rx_hat ( : , 2 ) =  Rx ( : , 2 ) . /sqrt (Rx ( : , 1 ) . *Rx ( : , 1 ) +Rx ( : , 2 ) . *Rx ( : , 2 ) ) ; 

%Signal  information 
c=3e8 ; 

W  tx=2*pi* (10e9) ;  %Carrier  Freq  is  lOGHz 

K  tx=W  tx/c; 

%Waveform  information 
fs  =  20e6; 

tl  =  0:1:20;  %pusle  transmit  time 

t2  =  0:1:1979;  %listening  time 

period  =  lOOe-6; 

T_period=0 : ts : lOOe-6; 

T=[tl  t2  tl  t2  tl  t2];  %period  of  lOOus 

y  =  rectpuls ( 0 , tl ) ; 
s  =  [y,  zeros ( 1 , length  (t2 ))] ; 
sp  =  [s  s  s]  ; 

S=fft(sp);  %FFT  of  signal 

w= (2*pi/period) *T; 

TT_Data= zeros ( 1 , length ( sp) ) ; 

%Generating  Target  Signal 

for  l=l:N_Tx  %For  all  Transmitter 

for  m=l:N  Rx  %For  all  Receiver 
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for  n=l:N  tt  %For  all  targets 

tau=  T_tx+ ( ( (Tx_mag (1,  : ) - 

(Tx_hat ( 1 ,  : ) *tt (n,  : )  '  ) )  +  (Rx_mag (m,  : ) - (Rx_hat (m,  :) *tt (n,  :)'))) /c) ; 
%time  delay 

phi=K  tx*Rx  mag (m, : ) - 

K  tx* (Tx  hat  (1,  : ) +Rx  hat (m,  : ) ) * (tt (n,  : )  +  ( (Rx  hat (m,  : ) * (Rx (m,  : ) - 
tt (n,  : )  ) ' ) *tt_vel (n,  : ) /c) ) ' ; 

alpha=l- (Tx_hat (1, : ) +Rx_hat (m, : ) ) * (tt_vel (n, : ) /c) ' ; 
TT_Data  =  TT_Data  + 

exp (i*phi) *exp (i*W  tx* alpha* T) . *ifft (S . *exp (-i*w*tau) ) ; 

end;  end;  end 
%Expected  Target  Position 

E  tt  x=-le2 : 5 : le2 ;  %Sampled  by  Range  Resolution 

E_tt_y=-le2 : 5 : le2 ;  %Sampled  by  Range  Resolution 

E_tt_y=f liplr (E_tt_y) ; 

%First  speed 

%Expected  Target  Velocity 

for  X  =  1:N  tt  %used  for  speed  calculation 

for  y  =  -200:200  %speed  of  target  in  x  direction 

for  z  =  -200:200  %speed  of  target  in  y  direction 

if  [tt_vel(x,l)  tt_vel(x,2)]  ==  [y  z] 

V  =  sqrt (y^2+z^2 ) ; 

if  V  >=  0  %function  to  find  speed  of 

if  V  <  50  %target  falls  in  range 

E_tt_vel_x=  y; 

E_tt_vel_y=  z; 

E_tt_vel_y=f liplr (E_tt_vel_y) ; 

E_tt_vel= [E_tt_vel_x  E_tt_vel_y] ; 

%% first  speed 

%Generating  Expected  Target  Database 
for  g=l : length (E_tt_y) 

for  h=l: length (E  tt  x) 

E_tt= [E_tt_x ( 1 , h)  E_tt_y ( 1 , g) ] ; 

E_TT_Data=zeros ( 1 , length (T) )  ; 

for  1=1 :N  Tx  %For  all  Transmitter 

for  m=l:N  Rx  %For  all  Receiver 

tau=  T_tx+ ( ( (Tx_mag (1,  : ) - 

(Tx_hat (1, : ) *E_tt (1, :)'))+ (Rx_mag (m, : ) - (Rx_hat (m, : ) *E_tt (1, :)'))) /c) ; 
%time  delay 

phi=K  tx*Rx  mag (m, : ) - 

K  tx* (Tx  hat(l,:)+Rx  hat (m, : ) ) * (E  tt(l,:)+((Rx  hat (m, : ) * (Rx (m, : ) - 
E~tt  (  1  ,  :y)  '  )  *E_tt_veI/c)  )  '  ; 

alpha=l- (Tx  hat(l,:)+Rx  hat (m, : ) ) * (E_tt  vel/c)'; 
E_TT_Data  =  E_TT_Data  + 

exp (i*phi) *exp (i*W  tx* alpha* T) . *ifft (S . *exp (-i*w*tau) ) ; 
end;  end 

I (g, h) =E_TT_Data*TT_Data'  ; 
end;  end 

for  n  =  1: length (E  tt  x)  %Thresholding 

for  m  =  1: length  (E  tt  y) 
if  C (n, m)  >=250 

D  (n,  m)  =  C  (n,  m)  ; 

end 

if  C (n, m)  <250 
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j=0; 

D(n,m)  =  j ; 

end;  end;  end 

A  =  max (max (abs (D) ) ) ;  %largest  value  in  D  used  for  normalization 
B  =  D./A; 

C  =  abs (B) ;  %Normalized  data 

for  n  =  1: length (E  tt  x) 

for  m  =  1: length (E  tt  y) 
if  D (n, m)  >=  0.9 

H(l,l)  =  E_tt_x(l,m); 

H(l,2)  =  E_tt_y(l,n)  ; 

H(l,3)  =  0; 
end;  end;  end 

stem3 (H ( 1 , 1 ) , H ( 1 , 2 ) , H ( 1 , 3 ) , ' . k' )  %image  of  data  created 

quiver 3  (H(l,l)  ,H(1,2)  ,H(1,3)  ,y,z,0,  'k'  ) 
axis([-100  100  -100  100  0  1]) 
grid 

end;  end;  end;  end;  end;  end 

S-5-2-5-S-5-2-5-S-5-5-2-2-S-5-2-2-S-5-2-2-S-5-2-2-S-5-2-5-S-5-2-2-S-5-2-5- 

ooooooooooooooooooooooooooooooooooooo 

%The  following  is  an  example  of  how  the  data  could  be  partitioned  for 
%multiple  velocities 

S-S-2-5-S-2-2-S-S-2-2-2-5-S-S-2-S-S-2-2-2-S-2-2-2-S-2-2-5-S-2-2-5-S-S-2-2- 

ooooooooooooooooooooooooooooooooooooo 

%Second  Speed 

%%Expected  Target  Velocity 

for  X  =  1:N  tt  %used  for  speed  calculation 

for  y  =  -200:200  %speed  of  target  in  x  direction 

for  z  =  -200:200  %speed  of  target  in  y  direction 

if  [tt_vel(x,l)  tt_vel(x,2)]  ==  [y  z] 

V  =  sqrt (y^2+z^2 ) ; 
if  V  >=  50 

if  V  <  100 

E_tt_vel_x=  y; 

E_tt_vel_y=  z; 

E_tt_vel_y=f liplr (E_tt_vel_y) ; 

E_tt_vel= [E_tt_vel_x  E_tt_vel_y] ; 

%Generating  Expected  Target  Database 
for  g=l : length (E_tt_y) 

for  h=l: length (E  tt  x) 

E_tt= [E_tt_x ( 1 , h)  E_tt_y ( 1 , g) ] ; 

E_TT_Data=zeros ( 1 ,  length (T) )  ; 

for  1=1 :N  Tx  %Eor  all  Transmitter 

for  m=l:N  Rx  %Eor  all  Receiver 

tau=  T_tx+ ( ( (Tx_mag (1,  : ) - 

(Tx_hat (1, : ) *E_tt (1, : ) ' ) ) + (Rx_mag (m, : ) - (Rx_hat (m, : ) *E_tt (1, :)'))) /c) ; 
%time  delay 

phi=K  tx*Rx  mag (m, : ) - 

K  tx* (Tx  hat(l,:)+Rx  hat (m,  : ) ) *  (E  tt(l,:)  +  ((Rx  hat (m,  : ) * (Rx (m,  : ) - 
E~tt  ( 1 ,  :T)  '  )  *E_tt_veI/c)  )  '  ; 

alpha=l- (Tx  hat(l,:)+Rx  hat (m, : ) ) * (E  tt  vel/c)'; 
E_TT_Data  =  E_TT_Data  + 

exp (i*phi) *exp (i*W  tx* alpha* T) . *ifft (S . *exp (-i*w*tau) ) ; 
end;  end 
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I (g,  h) =E_TT_Data*TT_Data' ; 
end;  end 

for  n  =  1: length (E  tt  x) 

for  m  =  1: length (E  tt_y) 
if  C(n,m)  >=  160 

D  (n,  m)  =  C  (n,  m)  ; 

end 

if  C (n, m)  <  160 
j=0; 

D(n,m)  =  j ; 
end;  end;  end 
A  =  max  (max  (abs  (D)  )  )  ; 

B  =  D./A; 

C  =  abs (B) ; 

for  n  =  1: length (E  tt  x) 

for  m  =  1: length (E  tt  y) 
if  D (n, m)  >=  0.9 

H(l,l)  =  E_tt_x(l,m); 

H(l,2)  =  E_tt_y(l,n)  ; 

H(l,3)  =  0; 
end;  end;  end 

Stems (H(l,l),H(l,2),H(l,3),'.k') 
quivers  (H(l,l)  ,H(1,2)  ,H(1,3)  ,y,z,0,  'k'  ) 
axis([-100  100  -100  100  0  1]) 
grid 

end;  end;  end;  end;  end; 


end 
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